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Abstract 

A model for S-wave rjir scattering is proposed which could be realistic in an 
energy range from threshold up to above one GeV, where inelasticity is dominated 
by the KK channel. The T-matrix, satisfying two-channel unitarity, is given in a 
form which matches the chiral expansion results at order exactly for the r/vr —)• 
rjTT, rjTT —>• KK amplitudes and approximately for KK —)• KK. It contains six 
phenomenological parameters. Asymptotic conditions are imposed which ensure a 
minimal solution of the Muskhelishvili-Omnes problem, thus allowing to compute 
the r]TT and KK form factor matrix elements of the 1 = 1 scalar enrrent from the 
T-matrix. The phenomenological parameters are determined such as to reproduce 
the experimental properties of the ao(980), ao(1450) resonances, as well as the 
chiral results of the t/tt and KK scalar radii which are predicted to be remarkably 
small at O(p^). This T-matrix model conld be used for a unified treatment of the 
rjTr final-state interaction problem in processes such as r/' —)• ijinr, (j) —)• or the 
T/TT initial-state interaction in r/ —)■ Svr. 
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1 Introduction 

The properties of the r]n scattering amplitude are much less known than those of nn 
or Kn scattering. In the framework of three-flavour chiral symmetry (in which the r] 
is a pseudo-Goldstone boson) a specihc prediction can be made that the r]n interaction 
should be considerably weaker than the nn or Kn interactions [1] at low energies. This 
feature has not yet been verihed either experimentally or in lattice QCD. It is possibly 
related to the apparent absence of a broad light 1 = 1 scalar resonance. 

A global description of nr] scattering (in particular of the elastic channel and the 
leading inelastic channel nr] —)■ KK) would enable one to perform a universal treatment 
of the hnal-sate (or initial-sate) interaction involving the nr] system. A particularly 
interesting application would be to the r] ^ 3n amplitude. Precision measurements 
of these decay modes should be exploited in an optimal way for the determination 
of isospin violating quark mass ratios. For this purpose, it is necessary to combine 
chiral expansion expressions with general dispersive treatments of rescattering [2, 3]. An 
extension of these approaches to include r]n rescattering would allow one to take into 
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account explicitly the Oq — /o “mixing” effect^, which was claimed to be signihcant [5] 
for T] —)■ Stt. 

The available experimental information on rjn scattering have been derived via the 
hnal state interaction effects in production processes and they concern, essentially, the 
properties of the resonances. The two prominent resonances which have been observed 
in the 5'-wave are the ao(980) and the ao(1450). We wish to address here the problem 
of determining more global properties of the S'-wave amplitude i.e. the determination 
of phase shifts and inelasticities in the small to medium energy range such as to be 
compatible with the properties of the resonances and also obey further theoretical con¬ 
straints. 

Several models of the r^vr S'-wave scattering amplitude have been proposed in the 
literature [6-9]. Our approach enforces a correct matching with the chiral expansion 
of the amplitudes at low energy in a way somewhat similar to refs. [6, 7]. In addition, 
we propose here to consider the form factor (and Fg^) associated with the scalar 
isovector current operator ud, in parallel with the T-matrix. Form factors are the 
simplest quantities to which analyticity based hnal-state interaction methods can be 
applied. We will follow the same general method which was proposed for the scalar 
isoscalar pion (and kaon) form factors [10] and proved capable of determining the scalar 
radius of the pion {r'^Yg rather accurately (see refs. [10-16] for theoretical calculations, 
and refs. [17-19] for lattice determinations). Its application to the strangeness changing 
Ktt scalar form factor and the corresponding scalar radius {r^g^ were discussed in 
refs. [20-22]. Form factors are constrained by chiral symmetry at low energy and, even 
though the convergence of the three flavour chiral expansion may be rather slow, one 
still expects correct order of magnitudes to be provided at order At this order, a 
simple relation between the pvr and the Kn scalar radii is predicted 


yyK 


0.52 ±0.02 . 


( 1 ) 


This relation implies that the pvr radius is remarkably small — 0-1 We will 

show that this result provides a stringent constraint in the determination of the phase 
shifts and inelasticities. 

The plan of the paper is as follows. We start with the chiral perturbation theory 
(ChPT) expansions of the scalar form factors Fg^, Fg^ and with the rjn and KK scat¬ 
tering amplitudes at next to leading order (NLO). Next, we recall the general dispersive 
integral equations from which one can compute the form factors starting from a given 
T-matrix, provided suitable asymptotic conditions are imposed. We then describe our 
chiral A'-matrix type model for the T-matrix, which involves six phenomenological pa¬ 
rameters. It is designed such that, at low energies, the contributions involving these 

^This effect was first discussed in ref. [4]. It can be seen as a superposition of the two physical 
resonances ao(980), /o(980) in the rjir —>■ tttt scattering amplitude. 
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parameters have chiral order p® (that is, NNLO) and that a proper matching with the 
ChPT expressions at NLO holds except, however, for the KK —)■ KK amplitude, for 
which the matching is only approximate. Finally, the determination of the phenomeno¬ 
logical parameters is discussed such as to satisfy the experimental constraints on the oq 
resonances and the chiral constraints on the scalar form factors. 


2 ChPT expansions of form factors and 

scattering amplitudes 

2.1 Form factors and scalar radii 

Let us introduce the following two form factors associated with the isospin one charged 
scalar operator ud 


Bq Fg 


BoF^^'is) = {ri{pi)Tr^{p2)\ud{0)\0) 

s) = {K\p,)K+{p2)\ud{0m 


KK. 


( 2 ) 


where s = {pi + P 2 )‘^- We have computed these form factors at next-to-leading order 
(NLO) in the chiral expansion. The detailed expressions are given in appendix A. From 
eqs. (52), (53) in that appendix, it is easy to derive the expressions of the scalar radii, 
which are defined as 

(3) 


= 6F|'«(0)/F|'«(0) 


PQl 


For pir and KK one obtains 
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(5) 


where Lp, Rpq are logarithmic functions of the pseudo-scalar meson masses. 


Lp = log 


m 


/i 


BpQ — 


w?p \og{w?p/rriq) 


nil 


m 


Q 


( 6 ) 


with p a renormalisation scale. These scalar radii depend on only one of the Gasser- 
Leutwyler coupling constants [23], L^. It is instructive to compare them with the anal¬ 
ogous Ktt scalar radius associated with the strangeness changing scalar current, which 
also depends only on Lg [24], 



6 

F3 


AU - ^ ^ 

5 8167r2 


QLk + ^Rttk + R 


i-nK 


(7) 
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The explicit expression of 62 , from ref. [24], is reproduced in appendix A. One remarks 
that the three scalar radii have exactly the same dependence on 

the coupling Lg, which means that they should be equal in the large Nc limit of QCD. 
In reality, they are rather different. Using e.g. Lg = (1.23 ± 0.06) • 10“^ (from ref. [25], 
see sec. 2.3 below) one hnds^ for r/vr and KK 


(r2)^^ = 0.092 ± 0.007 fm% 

(r2) = 0.136 ± 0.007 W , 

while for TCvr, one hnds, 

(r^) = 0.177 ± 0.007 W . (9) 

This shows that the rjir scalar radius is suppressed by a factor of two as compared to 
the Kti scalar radius. 


2.2 Scattering amplitudes at 0{p^) 


We consider the three scattering amplitudes involving the and the channels 

and we label the channel as 1 and the channel as 2. At chiral order the 

amplitudes read. 


Tg{s,t,u) 

Tg{s,t,u) 


9 

mt 

TT 

^ (3s - Ami) 

TT 

-^{s + it-u)) . 

TT 


( 10 ) 


The corrections of chiral order to these amplitudes can be expressed in terms of a set 
of functions of one variable, analytic with a right-hand cut, according to the so-called 
reconstruction theorem [26] (see also the review [27]), 


7(i;(s,f,n) =Ul\s) + Ul\u) + W^\t) 

T^A)i.s,t,u) = U^^{s)+ [W^^{t) + {s-u)Wi{t) + {t u)] (11) 

7(^J(s, t, u) = U^^(s) + (t- u)Ui(s) + Vo(t) + (s- u)Vi{t) + W^‘^{u) . 

The detailed expressions of the functions Uj, V) are given in appendix B. 

The resulting amplitudes are equivalent to previous calculations [1, 7]. We dehne the 
partial-wave amplitudes as 

Tf{s) = r (12) 

^The following input numerical values are used throughout this paper (all in GeV): m,r = 0.139568, 
mx = 0.4957, = 0.547853, Ur = 0.09221. 
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( 13 ) 


such that the unitarity relation, in matrix form, reads 


Im Tj{s) = Tj{s) S(s) Tjis) = T}{s) S(s) Tj{s 


with 


and 


ai(s) = 




(T2{s)6{s - Amj^) 


(14) 




(^2(8 = 


s — 4m| 


*' V(«) = + • 

(15) 


The relation between the partial wave S- and T-matrices then reads 

5j(s) = 1 + 2i^/^Tj{s) . 


(16) 


t, u{z^^) = i ( 2ml + 2ml - s± 


A2 

r)77 


In eq. (12), designate the cosines of the centre-of-mass scattering angles, which are 
related to the Mandelstam variables by 

•^T77r('5)'^ 

t, u{z^^) = I + ml + 2m\ - s± a 2 {s) z^'^^ 

t,u{z'^‘^) = |(4m|- - s)(l =F z^'^) 

with = rri^ — rri^. The hrst two of these relations become singular when s ^ 0. 
This implies that the chiral expansions of the 7]tt —)■ 7]tt and r]TT ^ KK partial-wave 
amplitudes become invalid when s is too close to zero. If we assume a domain of validity 
for the expansion of the unprojected amplitudes when |s|, |f|, \u\ < 0.5 GeV^, then the 
chiral expansions of the partial-wave amplitudes Tj^, should converge with s lying 
in the range 0.17 < s < 0.5 GeV^ and 0.05 < s < 0.5 GeV^ respectively. 

From now on, we will consider only the J = 0 partial-wave and will drop the J 
subscript. With the subscript now indicating the chiral order, the J = 0 partial-wave 
amplitudes at 0{p^) are simply derived from (10) 


T(l(s) = 


m: 


IGtt 3FJ ’ 


mis) = 


1 ^6(3, 


Am]^) 


IGvr 


12F2 


mm = 


IGtt AF^ 


(18) 


The corrections of chiral order to these J = 0 partial-wave amplitudes can be written 
as 


AF) = 


IGvr 


s) + uiiis] 


( 19 ) 


where 


U, 


11, 


s = 


'-1 


dz" (c/„“(«) + ir„”(i)) 
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( 20 ) 


= j (W'o’W + (s - m) 

U?M = 5 {V„(t) + (s-u) Vdt) + WS^(u)) 

The functions Ul\s) carry the left-hand cuts of the partial-wave amplitudes T*h These 
cuts are as follows [28]: 

A real cut on [—oo, (rUr, — rriTrY] a complex circular cut centred at s = 0 with 
radius 

A real cut on [—cxo, 0] and a complex quasi-circular cut which intersects the real 
axis at —Ar^T^rriK/{rrix + m,,) and Ar^T^niK/ijriK + 'm-,,). 

A real cut on [—cxo,4m|- — 4m^]. 

As a hnal remark, at NLO, each one of the functions Uq , Ui, Vj can be written 
as the sum of a polynomial part and one involving a combination of functions Jpq (see 
appendix B). The latter part is constrained by unitarity. For instance, for the functions 
one can write, in matrix form, 

U„(s) = P„(s) + (167r)2T,2,(s) T|2 ,(s) , (21) 

2.3 Influence of the 1/Nc suppressed couplings 



10^ L[ 

10^ 

10^ LI 

10^ 

10=^ 

10=^ 

10^ L^ 

0 

CO 

00 

(A) 

1.11 

1.05 

-3.82 

1.87 

1.22 

1.46 

-0.39 

0.65 

(B) 

1.00 

1.48 

-3.82 

0.30 

1.23 

0.14 

-0.27 

0.55 


Table 1: Two sets of central values of T[(p) with p = 0.77 GeV from NLO fits performed 
ref. [25]. 

The values of the low-energy couplings (LEC’s) L[, i = 1 • • - 8 are needed in order 
to evaluate numerically the chiral amplitudes. A recent update of the values of the 
couplings has been presented in ref. [25] based on global hts involving a number of 
low energy observables. We reproduce in table 1 two sets of values which correspond to 
NLO expansions (which seem appropriate here since we are using NLO formulae). The 
set labelled (A) in table 1 corresponds to an unconstrained ht and it leads to rather 
large values of the couplings L 4 , Lq and L 2 — 2Li which are suppressed in the large Nc 
limit [23]. The set (B) in the table corresponds to a ht which is constrained to enforce 
compatibility with the results from lattice QCD simulations on L 4 and Lg. We will 
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consider it to be more plausible, since the strong deviations from the large limit are 
possibly an artifact of attempting to reproduce certain observables which are sensitive 
to NNLO rescattering effects (like the J = J = 0 tttt scattering length) using NLO 
formulae. Fig. 1 illustrates the sensitivity of the J = 1 amplitudes considered here to 
the l/Nc suppressed couplings. The shape of the rjir —)■ rjTi amplitude is quite different 
if one uses the set (A) or the set (B). This is also reflected in the values of the J = 0 
threshold parameters. Dehning the scattering length oq and the scattering range Bq as 
in ref. [ 1 ], 

—^ T^^(s) = Oo + &o + • ■ • (22) 

Vs 

with ^/s = ml + one hnds 


tUtt oo = 6.7 ■ 10 ttIt, 6o = —15.0 ■ 10 ^ (Large L 4 , Lg) 

= 16.2 ■ 10“^, 60 = 10.6 • 10“^ (Small L 4 , Lg) . 


(23) 


The two sets of couplings thus lead to rather different values of the scattering length oq 
while the values of the scattering range 60 differ in their sign. At leading chiral order, 
one has oq = 6.2 ■ 10“^, 60 = 0. At NLO, a low-energy theorem (LET) for og was 
derived in ref. [29], in the form of a linear relation 

~ ®0,7r7r|7V2.O ^ 

where Oq is the tttt scattering length with J = 0, / = 2 and A, /i are simple functions 
of the masses and the decay constants Ft,, Fk- The most precise deter¬ 

minations of the S'-wave tttt scattering lengths are based on Roy equations solutions. 
Using the values quoted in two recent analysis of these equations [14, 30] in the LET 
relation (24), one obtains 

=-0.0444 ± 0.0010 (ref. [14]) —^ oq = (- 0.22 ± 6.26) ■ lO-^ 

= -0.042 ± 0.0040 (ref. [30]) —> ao = (14.8 ± 25.0) • 10 -^ . 

This illustrates that the LET is practically useful only if Oq t,,, is known to a very high 
precision. The result of ref. [14] is associated with a rather small error of 2.5%. However, 
the result derived from the Roy equations concerns the physical value of the scattering 
length rather than the NLO value which enters into the LET. An additional error should 
therefore be introduced in eq. (25) in order to account for the difference Oq — Oq 
which could easily be as large than 5%. This observation then limits the effectiveness 
of the LET for determining og. 

The partial-wave amplitude vanishes at s = 0 at leading chiral 

order (18). This zero, however, is accidental since it is not associated with a soft pion 
theorem. Fig. 1 shows that, indeed, the NLO corrections are substantial. The corrections 
corresponding to the Li set (B), with small 1/Nc violations, have a more reasonable size 





s [GeV^] 


Figure 1: Real parts of the three J = 0 partial-wave amplitudes r/vr'*' —>• r]7r~^, T/vr"*" —)• K^K~^ 
and K^K~^ —)■ at leading and next-to-leading order in ChPT. 


than those from set (A). The amplitude r]Ti~^ —)■ has a zero at s = 4m|-/3 at 0{p^) 

which corresponds to a soft pion Adler zero. Fig. 1 shows that the NLO corrections are 
rather small in this case and that there is little difference between the couplings of set 
(A) and set (B). 

3 Form factors from dispersive integral equations 

We follow here a general approach to the construction of form factors which implements 
unitarity relations and chiral constraints and, additionally, impose the absence of zeros 
and consistency with the QCD asymptotic behaviour. We will briefly review this method 
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below, which was applied previously to the scalar tttt and tiK form factors [10, 20], and 
allows one to relate the form factors and the corresponding S'-wave scattering amplitudes 
via a set of integral equations. The J = 1 scalar form factors , Fg^ which we will 
discuss here were considered previously in ref. [31]. The approach followed in ref. [31] 
differs from ours in that the constraints on the zeros and the asymptotic behaviour were 
not imposed. 


3.1 Phase dispersive representation 

The crucial property of two-meson form factors is that they can be dehned as analytic 
functions in the complex energy plane, with a cut lying on the positive real axis in the 
range s > {mp + mq)'^ [32]. In the asymptotic region, |s| oo, the general arguments 
concerning exclusive processes in QCD [33] predict that a two-meson scalar form factor 
Fs should obey a power law behaviour. 


Cs(s)L^oo ~ 1/s 


(26) 


up to logarithms. Making the assumption that the form factor Fs has no zeros in the 
complex plane, one can derive a minimal phase dispersive representation (e.g. [34]), 


Fs{s) = Fs{0) exp 


s 

TT 


r Ms') 

Jso s'{s'-s) 



(27) 


where the phase is dehned from F{s -|- ie) = |T)j(s)| exp{i(j)s{s)). The QCD asymptotic 
behaviour (26) is reproduced from eq. (27) provided that the phase has the asymptotic 
limit: 

lim (j)s{s') = IT . (28) 

s'—>-+C» 

The scalar radius, hnally, is given by a simple integral as a function of (ps, 



6 n Ms') 
TT Jso (s'y 


ds' . 


(29) 


If n complex zeros were present, then the right-hand side of eq. (27) would have to 
be multiplied by a polynomial of degree n and the asymptotic phase would have to be 
(n -|- l)7r. The minimality assumption is equivalent to stating that the increase of the 
phase in the energy region i/s > 2 GeV should be less than vr. This is plausible since 
no sharp resonances are present in this region. 


3.2 Determination of the form factors from the T-matrix 

As emphasised in ref. [15], these phase relations are of particular interest for those form 
factors which involve at least one pion, Fg^ with P = n, K ot r] which interests us here. 


10 






This is simply because the scattering amplitudes ttP —)■ ttP are elastic in a finite low 
energy region. In this region, the form factor phase is constrained from Watson’s 
theorem to be exactly equal to the elastic scattering phase shift. The energy region 
in which inelasticity can be neglected to a good approximation extends up to the KK 
threshold for tttt and we expect the same property to hold also^ for irr]. The asymptotic 
value of the form factor phase is also known and one may estimate that should be 
smoothly approaching its asymptotic value when ^/s ^ 2 GeV. There only remains to 
determine in the intermediate energy region that is, in the case of r^vr, in the region 
1 < ^/s ^ 2 GeV. In this region, we further expect that the fastest energy variation 
should take place close to 1 GeV, associated with the sharp onset of inelasticity triggered 
by the presence of the ao(980) resonance which is known to couple strongly to the KK 
channel [35]. This suggests to consider a framework which takes into account only the 
dominant inelastic channel and ignores all the other ones. In this case, the two form 
factors Fg^, Fg^ obey a closed set of Muskhelishvili-Omnes coupled integral equations, 

fFr{s)\_l r ds' fT^\s) T^^is)y f F^g^s') \ 

KFg^Hs)) tt V - 5 T^\s)) F^^{s')e{s' - Ami)) ' 

(30) 

These equations encode the property of analyticity of the form factors, the asymptotic 
behaviour (which allows for an unsubtracted dispersive representation) and two-channel 
unitarity. One can express the two-channel S'-matrix in terms of two phase shifts and 
one inelasticity parameter in the usual way. 


^ %\/\ — "tf 

We assume the following asymptotic conditions on the S'-matrix parameters 

lim T]{s) = 1, lim (5ii(s) S 22 {s) = 27r , 

S^OO S^OO 


0<r]<l 


(31) 


(32) 


which ensure that the so called Noether index [36] (see also [37]) associated with the set 
of singular integral equations (30) is equal to two. This, in general, implies that a unique 
solution is obtained once two arbitrary conditions are specified, for instance the values 
at s = 0: Fg^{0), Fg^{0), and that the solution form factors behave asymptotically as 

l/s |37]. 

In summary, solving the set of eqs. (30) for the form factors Fg'", Fg-^, one obtains a 
phase (pl which correctly matches with both the low and high energy limits expectations 
and provides an interpolating model in the intermediate energy region. The phase (pg^ 

^The inelastic mode -qw —>■ Stt is allowed already at threshold but the S-wave projection vanishes 
by parity conservation (since = 0“ for the Stt state). The modes qir —>■ Stt, qn —> qSir are strongly 
suppressed by phase space below one GeV. 
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is also provided. In this case, however, there is no constraint from Watson’s theorem at 
low energy. One expects that the form factor will be more sensitive than to 
the influence of the neglected inelastic channels. 

More generally, one can use the system of equations (30) to dehne the Omnes matrix 
which generalises the usual Omnes function [38]. Such a generalisation was hrst 
discussed in the case of tttt — KK scattering in refs. [39, 40]. The hrst column of 
the Omnes matrix is obtained by solving the system with the boundary conditions 
O^^(O) = 1, r2^^(0) = 0 and the second column by solving with the conditions r2^^(0) = 0, 
r2^^(0) = 1 (see in ref. [13] an appropriate numerical method for solving the linear 
system). The Omnes matrix allows one to treat the hnal-state interaction problem 
taking into account inelastic rescattering. For instance, one can express the 1 = 1 scalar 
form factors in terms of the 17 matrix. 




(33) 


4 Two-channel unitary T-matrix parametrisation with 
chiral matching 

We seek a parametrisation of the J = 0 T-matrix which; a) should satisfy exact elastic 
unitarity below the KK threshold and exact two-channel unitarity above, b) should 
correctly match with ChPT for small values of s , i.e. 

T«(s) - (T|(s) + (s)) = 0(p'=) . (34) 

and c) should be reasonably simple and flexible and be able to describe scattering in 
the low to medium energy region up to, say >/i ~ 2 GeV. We choose a representation 
somewhat similar to that proposed in ref. [41] to describe J = 0 nK scattering, belonging 
to the family of “unitary chiral” approaches. Such approaches were proposed, in the 
context of ChPT, hrstly in refs. [42, 43] and multichannel extensions were discussed in 
refs. [44, 45] (we refer to the review [46] for a survey and a complete list of references). 
There are, however, some drawbacks to these methods. Poles can occur on physical 
sheets and, furthermore, the structure of the left-hand cuts is not quite correct. In 
particular, the left-hand cut of the chiral KK —)■ KK amplitude T^ 4 ^(s), which extends 
up to s = 4(m|- —m^) is propagated to the amplitude via the unitarisation method, 
which actually spoils the unitarity of in the elastic region. While the resulting 
unitarity violation is numerically small [7, 47], we will prefer here to maintain exact 
unitarity at the price of relaxing the matching condition for the component 
We start from a iP-matrix type representation for the two-channel T-matrix 

. ( 35 ) 
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This form is compatible with the symmetry of the T-matrix (*T = T) provided both K 
and $ are symmetric matrices. The matrix $(s) must also satisfy 


Im [$(s)] 


/ 9{s - {nir^ + m^)^)ai{s) 

V 0 9{s 


0 ^ 

4ml)a2{s)J 


(36) 


which ensures that the T-matrix obeys the unitarity condition, provided that the matrix 
K{s) remains real in the range (m^ -|-m 7 r)^ < s < oo. We take a representation of $(s), 
satisfying eq. (36), which is diagonal and contains four phenomenological parameters 


^ + 167rJ^7,(s) 


0 






(37) 


The parameters Oj, /9j are assumed to be 0(1) in the chiral counting. The ii'-matrix is 
written in terms of components with a dehnite chiral order. 


K{s) = K( 2 )is) + K( 4 )(s) + K(6)(s) (38) 

where, as before, the subscript denotes the chiral order. In order to satisfy the matching 
condition (34) one must have, 

K(2 ){s) = T(^2){s), T(4)(s) = 1T(4)(s) T(2)(s)$(0)(s)T(2)(s) . (39) 

One can then express ^T( 4 ) in terms of the polynomial and left-cut functions dehned 
from eqs. (19) (20) (21) (see also appendix B) 

*:(4)(i>) = ^ (p„(s) + c/„(s)) - rp,(ii) T(2 i(s) . (40) 


As explained above, we must use an approximation to the function 1 /q^ which has no 
cut on the real axis in the range s > (m^ -|- m^)^. This may be done by removing the 
parts which are proportional Jmrit) and Jrinit) (see eqs. (74)) from the two functions 
Vo(t) and Vi(^), which appear in the angular integral which gives (see eq. (20)). 
Figure 2 compares this approximation of to the exact function. 

Finally, lT(6)(s) is taken to be a pole term with the part removed. 



9i9j f 1 _ 

IOtt — s ml 


(41) 


We model the couplings gi, g 2 such that they behave as based on a scalar reso¬ 

nance chiral Lagrangian analogous to the one introduced in ref. [48] 

9i ml) + 2c'^ml) , 
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S (GeV)2 

Figure 2: Comparison of the real part of function ?7 q^ with the approximation used in the 
unitary representation (35). Also shown is the imaginary part of 

92 =-4(Cd (s - 2m^) + 2c'^ m^) . (42) 

■^TT 

We will discuss in sec. 5 how the phenomenological parameters may be determined from 
experimental information on the properties of the ao(980), ao(1450) resonances as well 
as chiral constraints on the amplitudes and on the 1 = 1 scalar form factor. Figure 3 
illustrates how the unitary amplitudes parametrised as described above correctly match 
with the NLO chiral amplitudes at low energy. 


5 Phenomenological determination of the phase shifts 
and inelasticity and the 1 = 1 scalar form factor 

5.1 Experimental information on irr] irr] and ttt] KK 
scattering 

Let us hrst consider the irri —)■ tit] amplitude below the KK threshold. In this region, 
rjTr scattering should be approximately elastic. The tit] scattering phase shift below 1 
GeV should be controlled by the values of the threshold parameters Oq, &o on fhe one 
hand and the properties of the ao(980) resonance on the other. We will consider that the 
values of oq, foo corresponding to the set of L's with small L 4 , Lq (set (B), see table 1) 
are the most plausible. In this case, oq and bo are both positive and one expects that 
the phase shift will be positive in the whole elastic region. A different possibility was 
investigated in ref. [49]. 
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Figure 3: Comparison of the real parts of unitary partial-wave amplitudes given from 
eq. (35) and the corresponding chiral amplitudes at NLO. 


The ao(980) is a well established resonance but its shape is not well described by 
a simple Breit-Wigner form because of the vicinity of the KK threshold. This partly 
explains the dispersion in the values of the mass and width quoted by the PDG [50]: 

= 980±20 MeV, Tag = [50 — 100] MeV. A comparison of a number of determinations 
of the Til amplitude near the KK threshold based, in particular, on the popular Flatte 
model [51] is performed in ref. [52]. The corresponding rjn phase shifts are plotted on 
Fig. 10 of that reference, from which one can deduce that the value of the phase shift 
at the KK threshold lies around 90°, 

5ii(2mi^) = (90±20)° . (43) 

This is also satished in the models of refs. [8] and [9] which give, respectively, 5ii{2mK) = 
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95° and (5ii(2mx) = 77°. 

The ao(980) resonance corresponds to poles of the amplitude in the complex plane on 
the second and on the third Riemann sheets which can both be near the physical region 
since the mass is very close to the KK threshold. For dehniteness, we will rely here on 
the recent determination by the KLOE collaboration [53]. It is based on measurements 
of the 0 —)■ decay amplitude with both high precision and high statistics. Based on 
the best £t performed in ref. [53] (using the theoretical model from ref. [54]) the location 
of the poles can be deduced to be 


^5^(980) = (994 ±2-z(25.4±5.0))MeV 

(958 ± 13(60.8 ± 11.5)) MeV. 


(44) 


In the [1 — 2] GeV energy region, a second resonance, the ao(1450), first reported in 
ref. [55] was later identified in pp decays at rest (e.g. [56-58], see also [59] who re¬ 
analysed the data). This resonance should correspond to a pole on the third Riemann 
sheet. Based on the value of the mass and width quoted in the PDG, we can estimate 


= (1474 ± 19 - 2 (133 ± 7)) MeV . (45) 

A further property of the ao(1450) is that it has approximately equal decay widths into 
717] and into KK. We will implement this feature by requiring that the J = 0 cross 
sections for pir —)■ pir and r/vr —)■ KK should be approximately equal when y/s = 1.474 
GeV. In our two-channel framework, these cross sections have the following expressions 
in terms of the phase shifts and the inelasticity parameter 

(T(r77r ^ t^tt) = -r- [77 — l|^ , aipn ^ KK) = (l — rf) (46) 

Prj'K trri'K 

and we expect that p should reach a minimum at the mass of the ao(1450) resonance. 
If the minimum is close to zero, the two cross sections will be approximately equaP. In 
this situation, we expect a rapid variation of the phase shifts (fn, 622 (possibly becoming 
discontinuous if r; = 0) at the energy y/s = 7nao(i45o). In contrast, the sum of the two 
phase shifts (which is also the phase of 512 ) should be a smoothly varying function. It 
is convenient to characterise the global behaviour of the A-matrix in the [1 — 2] GeV 
region in terms of the value of this phase sum (fn -|- 612 when y/s = mag( 1450 ) 

<5i 2 = <5ii(x/i) + . (47) 

Let us now return to the parametrisation of the T-matrix described in sec. 4. The 
T-matrix elements in this model have analyticity properties and can be dehned away 

^Equality of the two cross sections occurs either when 77 = 0 or 77 = cos 2 ^ii 
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Figure 4: Phases hn, (522) their sum and the inelasticity ij from the T-matrix model of sec. 4 
corresponding to several imposed values of 812 (defined in eq. (47)). 


from the physical region, in the complex energy plane. Using eq. (35), the poles of the 
T-matrix correspond to the zeros of the determinant 

A(s) = det[l - K{s)^{s)] . (48) 

Recalling that the extension of the loop functions Jpq to the second Riemann sheet are 
dehned as 

J’pqM = Jpq(.s) + ‘ (49) 

then, the extension of the T-matrix elements to the second Riemann sheet is performed 
by replacing Jr^n{,s) by in the matrix $. Similarly, the extension to the third 

Riemann sheet is performed by replacing both and by and in $. 
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This T-matrix model involves the phenomenological parameters: ai, /9i, / 52 , mg, 
c'^, c'^. For simplicity, we will keep the ratio c^/c'^ hxed and allow only six parameters 
to vary. We determine them by imposing six conditions on the T-matrix: 

a) As hrst four conditions, we impose that the real and imaginary parts of the poles 
SaJ( 980 ) and s"fi 45 o) be reproduced. 

b) As a hfth condition, we impose that the minimum of the inelasticity parameter at 

= rnao(i 45 o) be close to zero (in practice, we used ?7mm ~ 0.05, as in ref. [60]). 

c) As a hnal condition, we choose a value for the phase 612 as dehned in eq. (47). 

Within this model, having imposed the hrst hve conditions, the value of 812 is found 
to be bounded from above: 612 ^ 205°. In addition, consistently with our assumption 
that most of the phase variations should take place below 2 GeV, it seems plausible 
that the phase sum (5ii -|- ^22 should not be smaller than its value at the mass of the 
ao(980), i.e. one should have 812 ^ 90°. Fig. 4 shows results from this model for the 
phases ^n, 822 and the inelasticity rj as a. function of energy, corresponding to several 
different imposed values of 812 - One observes that the two phases ^n, 822 undergo a 
sharp variation, in opposite directions, close to the mass of the ao(1450) resonance. The 
hgure illustrates a pattern where (5ii increases while ^22 decreases. However, a small 
modihcation of the phenomenological parameters which enter into the T-matrix model 
can lead to a pattern with a reversed behaviour (with (5ii decreasing and ^22 increasing) 
which would then be similar to the one obtained in ref. [60]. In contrast, the phase sum, 
811+822 is completely stable and always increases smoothly as an effect of the resonance. 
This ambiguity, which can be viewed as a ±7r ambiguity in the individual definition of 
5ii and 822 does also not affect observables, in particular, the determination of the form 
factors. 


<^12 

0.1 

02 

(3i (GeV-^) 

(82 (GeV-^) 

mg (GeV) 

A 

200° 

0.6265 

0.0988 

0.2495 

0.1476 

1.0571 

0.5704 

175° 

0.7427 

0.0781 

0.3085 

-0.0590 

1.0913 

0.8176 

150° 

0.8444 

0.0467 

0.2773 

-0.2085 

1.1258 

1.1017 

125° 

0.8765 

0.0016 

0.2134 

-0.3606 

1.1834 

1.6856 

100° 

1.0993 

-0.5055 

-0.0358 

-0.2722 

1.5130 

5.7024 


Table 2: Parameters of the T-matrix model corresponding to five fixed conditions (see text) 
and several input values of the phase 812 ■ The parameters c(„, c'^ are given in terms of A by 
(^m = Ac[J/ 2 with c[J = 28 MeV. 

Numerical values for the set of six parameters a,, /?,, mg, c'^ corresponding to several 
input values of 812 in the range 90° < 812 < 205° are given in table 2. The T-matrix is 
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v/s (GeV) 


Figure 5: Cross sections for rjir —)• t/tt and jyvr —)• KK in the vicinity of the ao(980) resonance 
from the T-matrix model, depending on the input value of 812 ■ The arrows show the integration 
limits used to define the branching fraction (50). 


not very sensitive to the value of the parameter c'^. Very similar results are obtained 
if one sets c'^ = 0 or = c'^. The numerical results shown in the table correspond to 
taking c'^ = c'^/2. In this model, the pole of the iC-matrix corresponds to two physical 
resonances. Table 2 shows that the mass parameter of the pole, ms, varies between 1 
and 1.5 GeV, while the value of the parameter c'^ varies in a rather large range from 16 
to 160 MeV, depending on the input value of the phase 612 - 

The properties of the ao(980) resonance (apart from the pole position on the second 
Riemann sheet which is held hxed) depend on the value of 612 - Figure 5 shows the 
two cross sections an-w^n-K, c^ti-k^kk iii fhe vicinity of the ao(980) resonance peak. We 
estimate the branching fraction in a simple way in terms of 

integrals over these cross sections 


BkK/tjtt 


r.E+ 


a. 


rf'K^KK 


{E) dE 



^?77r—>-777r (-S') dE 


(50) 


with E^ = mao ^ao- In this formula, we set ma^ = 988 MeV, which corresponds 
to the resonance peak in the cross sections and Fq^ = 50.8 MeV corresponding to 
twice the imaginary part of the pole position. We collect in table 3 the results for the 
branching fraction corresponding to different input values of 812 - The agreement with 
the experimental average quoted in the PDG, = 0.183 ± 0.024 is qualitatively 

reasonable, in particular for the smaller values of 812 - We also indicate in the table the 
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positions of the ao(980) pole on the third Riemann sheet (recall that the pole position 
on the second Riemann sheet is hxed), which is seen to move away from the real axis as 
612 is decreased. 


^12 

BriTT/KK 

v/iJF(MeV) 

200° 

0.095 

1022 62 

175° 

0.127 

1020 93 

150° 

0.148 

1009 - i 129 

125° 

0.170 

972 - i 192 

100° 

0.187 

749 - i 376 


Table 3: Some properties of the ao(980): values of the rjii/KK branching fraction and position 
of the pole on the third Riemann sheet depending on the input value of the phase 612 - 


5.2 Scalar form factors and the tjtz scalar radius 

In order to solve the integral equations (30) we must also dehne 5ii(s), 522('S), ri{s) 
for energies above the mass of the ao(1450) resonance such that the asymptotic condi¬ 
tions (32) are satished. For this purpose, we dehne a mapping u{s) such that 0 < u < 1 
when Si < s < 00 and then perform simple polynomial interpolations of the functions 
<^11, <^22, h ill terms of the variable u (see appendix C for more details, in practice we 
used =1.8 GeV). For a given value of the phase ^12, the T-matrix is completely 
specihed and one can derive the two scalar form factors by solving eqs. (30). 

The form factors turn out to be rather sensitive to the value of 812 . Fig. 6 illustrates 
the numerical results for the phase of the rjir scalar form factor, corresponding 
to different input values of ^12. The phase of the form factor displays a dip located 
in between the two Oq resonances. This behaviour is qualitatively similar to the one 
observed for the scalar form factor phases in the cases of the tttt or Ktt. A detailed 
discussion can be found in ref. [61]. The phase (fPg displays a bump, before the dip, which 
disappears when the input value of 812 is smaller than ~ 130°. Given the phase integral 
representation (29), we expect the r] 7 i scalar radius to decrease when 812 decreases. 
Numerical values of the scalar radii for the r^vr and the KK form factors are displayed 
in table 4 for given values of 812 in the range [100° — 200°]. In all cases, the dispersive 
result for exceeds the 0{p^) chiral value (8) (the same also holds for the KK scalar 

radius). However, one must also take into account the chiral corrections of order p® (or 
higher), the typical size of which can be as large as 20—30%. In the dispersive evaluation, 
even if the T-matrix elements were known exactly below 2 GeV, an error would arise 
from the asymptotic region. This is easily seen from the phase integral expression (29). 
The contribution to the pvr scalar radius from the integration region region > 2 GeV 
is relatively large ~ 30% and this could generate an overall uncertainty for {r‘^)Y of tho 
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Figure 6: Phase of the form factor obtained from solving the integral equations (30) with 
several input values of the phase 812 (see eq. (47)) in the T-matrix. 


order of 15%. The conclusion, then, is that the chiral result and the dispersive evaluation 
can be perfectly compatible provided the phase 812 lies in the following restricted range: 
90° < (5i2 < 125°. 


<^12 

200° 

175° 

150° 

125° 

0 

0 

0 

(fm^) 

0.185 

0.176 

0.166 

0.150 

0.122 


0.253 

0.248 

0.245 

0.233 

0.209 


Table 4: Results for the scalar radii obtained from solving eqs. (30) for the form factors 
depending on the input value for the phase 812 - 

Finally, £g. 7 shows the absolute values of the form factors ■ The size of 

the peak associated with ao(980) resonance is seen to be sensitive to value of the phase 
612 . We have verihed that the associated spectral function agrees with the one given in 
ref. [31] in the energy range s < 1.5 GeV^ when 812 — 100°. 


6 Conclusions 

We have proposed a model for the rjTr scattering T-matrix in the F-wave which satishes 
elastic unitarity below the KK threshold and two-channel unitarity above. The model 
is constrained by experimental inputs on the properties of the two resonances ao(980), 
ao(1450) and by chiral symmetry at low energy. In the simple JF-matrix type framework 
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Figure 7: Absolute values of the form factors (left) and (right) computed from our 
T-matrix model, corresponding to two input values of the phase 612 - 


which we have adopted it is possible to match correctly the two amplitudes titt rj-ir, 
rjTr —)■ KK with the chiral expansion at NLO while in the case of KK —)■ KK, the 
matching is only approximate (see sec. 4). Such a T-matrix could be realistic in an 
energy range ^/s ^ 1.3 GeV, where the inelasticity is effectively dominated by the KK 
channel. Formally, however, it is convenient to extend the model up to infinite energies 
such as to allow for a minimal solution of the associated Muskhelishvili-Omnes problem. 

A specific prediction of three-flavour ChPT is that the J = 0 rjir ^ titt scattering 
length is very small while the scattering range vanishes at leading order. The detailed 
predictions for these quantities at NLO are very sensitive to the values of the couplings 
T4, Lq which are l/Nc suppressed. We have used here the values of T4, Lq which are 
favoured by lattice QCD simulations. It would be a particularly interesting test of the 
chiral expansion, obviously, to have a verification of the t]ti scattering length also from 
lattice QCD. 

A supplementary chiral constraint which we have used is associated with the t]ti scalar 
isovector form factor. We have computed this scalar form factor from our two-channel 
T-matrix by solving the relevant Muskhelishvili-Omnes integral equations. While this 
model ignores other relevant inelastic channels (like ttt]') it is nevertheless plausible that 
it should be able to describe how the phase of the form factor behaves in approximately 
the same energy range where the T-matrix is realistic. Above this point, the model 
simply serves to interpolate the form factor phase monotonically towards its known 
asymptotic value. We find that the small value of the rjTr scalar radius in ChPT at NLO 
can be understood in this approach and that this requirement constrains the increase of 
the sum of S'-matrix phases (5ii -|- 622 in the 1 — 2 GeV energy region. One should keep 
in mind the uncertainties on the size of the NNLO effects on the ChPT side and those 
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from the energy range above 2 GeV on the dispersive side. The dispersive calculation 
suggest that the NNLO corrections to {r‘^)Y should tend to increase its size. It would 
again be extremely useful to have results from lattice QCD for this quantity. 

The computation of the 2x2 Omnes matrix 17 is a straightforward extension of the 
form factor calculation. In principle, the 17 matrix allows one to treat the rjir rescat¬ 
tering effects in a unihed way, in a number of processes for which recent measurements 
have been performed like rj' —)■ r^Trvr, 0 or 77 7 ]tt. The consideration of rjTr 

rescattering is also necessary in the case of the 7 —)■ Stt amplitude in order to account 
for oo — /o mixing within a dispersive approach. The rjTr scalar form factor itself appears 
in the isospin suppressed r —)■ r]7ru amplitude, along with an electromagnetic induced 
scalar form factor (and the vector form factor). This decay mode has not yet been 
observed but could possibly be studied at the super-B or future charm-tan factory. 
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A The 1 = 1 scalar form factors at NLO 


We consider the two form factors defined in eqs. (2). At leading order of the chiral 
expansion, the form factors are simply constant, 


fs’m = -Fs™(0) = 1 (iO) 


(51) 


Computing and adding the next-to-leading order corrections, the form factors can be 
written as 


FTM = -Fr(O) p + ^ 




4Lr + 

h I ^ 

TT 

1 ,. n „ . ^ . m. 


4F2 


(4m^ - 3 s) Jkk(s) + 


^p2 '^Fri\S) 


(52) 


^ 4L5 -h ~ 4) + 2 Lr^ + Lk) 


Jkk{s) 


4F2 


m 


{4:712^ - 3s)J^^(s) > , 


(53) 


where we have introduced the notation 


m 


Lp — log — p, RpQ — 


mp log(mp/''^Q) 


/i 


mp — m 


(54) 


Q 


and Jpq{s) are the loop functions dehned to vanish at s = 0 (we use the same notation 
as ref. [23]), 


■Jpq{s) — 


ds 


, V^pq(^') 


167r2 j^rnp+mg)'^ (s')^(s' - s) 


IhTT^ - 


1 + 


jPQ 

^PQ 


^pq\ , mp \/^pq{s) SpQ — s + \/\pq{s) 


log — + 

mq 


2 s 


log 


SpQ — S — \/Xpq{s)- 


(55) 


with 


FpQ — trip + ttiQ , ApQ — tpp — ttiQ , \pq{s) - - 2EpQ s + ApQ 


( 56 ) 


The expression for Fg^(O) is given by 


FFiO) 0 { 


1 + 


m 


K 


F3 L 


32 

- 64 + 32 L;: - — - 16 LI 

o 
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^(1-?L, + 2L^) 

16 

32 + 64 + 16 - 8 LI 


m: 

F2 L 


^7 


^6 


1677^ 


1/? -^r -^r 

■3 18 2 " 


(57) 


and the expression of Fg(0) reads 


2 r 


Ff *-(0) =1 + ^ [16 (2 LS - L3 + 16 (2 LJ - L;) + ^ (I fl., + “ L,) 


^ P '2 




9 


(58) 


A.l Remarks on Fg^{0) 

The value of Fg^{0) can be simply related to the mass difference. Indeed, 

using isospin symmetry, on can express the form factor Fg^ as 


BoF^^is) = 


,K^K^ - K+K- ,uu-dd 


V 2 


V 2 


| 0 ) 


(59) 


Then, writing the quark masses as 


= m -A, 


dU) fh T 


Feynman-Hellman’s theorem yields the following relation, 

d 


BoF#'^(0) = — (Afi. - Afjq . 


du 


(60) 


(61) 


One can easily reproduce eq. (58) using this relation and the chiral formula for the 
mass difference M‘^o — from ref. [23]. Using this formula, one can also derive an 
alternative expression for ^^-^(0), 


F#*(0) 


(m\-ml) r2 + 1 

/ /v \ 7-» ^ 

[mg — m)Bo r + 1 


(62) 


where r is the quark mass ratio mg/m and = 2m\^/m^ — 1 is the value of this ratio 
at chiral order p^. The deviation of the value of F|^^(0) from 1 can thus be interpreted 
as a measure of the size of the 0 {p^) corrections in the chiral expansion of the mass 
difference m\ — m^. Table 5 below shows that, if one uses the set of L[s with large L 4 , 
Lg, this correction is rather large (of the order of 40%). 
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We can also perform a verification of the value of Fg^[0). Using the Ward identity 
in pure QCD, 

id^u'j^d = {md — mu) ud (63) 

we can relate to the value at zero of the rjTr vector form factor (normalized 

as in ref. [62]) when = 0 




V^(m2 - ml) 
{md - mu)B q 



(64) 


Inserting the chiral expansion expressions for ml from ref. [23] and /+^(0) from 
ref. [62] one can recover eq. (57). 

The numerical values of Fg^{0), F^^{0) are needed as input in order to solve the in¬ 
tegral equations (30) for the scalar form factors. The values at s = 0 are rather sensitive 
to the l/Nc suppressed couplings L 4 , Lq as can be seen from table 5 below. However, 
the determination of the scalar radii from the integral equations depends 

only on the ratio Fg^{0)/F^^{0). It is easy to verify that this ratio is independent from 
L 4 , Lq at NLO. 



0 {v^) 

Small L 4 , Lg 

Large 

fy(o) 

0.816 

0.826 

1.421 

F#^( 0 ) 

1 

0.816 

1.428 


Table 5: Numerical values of Fg^{0), Fg^{0) in the chiral expansion at LO and at NLO using 
two sets of low-energy couplings (see table 1). 


A .2 Expression of 62 


We reproduce here the detailed expression (as given in eq. 6.2 of ref. [24]) for the term 
62 which appears in the chiral expansion of the Ktt scalar radius at order (see eq. (7)) 


with 


52 = 


1927r2F2 L 


h2{x) = 


15 


mz 


m 


K 


+ 


19 mj^ + 3m 


mj^ -|- m' 




mZ 


m 


K 


3 /1 -|- X 


2 V 1 


X 


3x (1 -b x) 

+ ^ log(x) . 


(1 


X 


18 


(65) 


( 66 ) 


B NLO contributions to / = 1 scattering ampli¬ 
tudes 


We give below the expressions of the chiral NLO contributions to the one-variable func¬ 
tions associated with the amplitudes pTr’*' —)■ r] 7 r~^, r] 7 r~^ —)■ and . 
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B.l The rjTT'^ amplitude 

The O(p^) part of the amplitude was written in terms of the two functions Uq^, 
(eq. (11)). They can be expressed as follows 


3 1 


4 m + Ly3) - (1 + Lk) 


8 IhTT^ 


+ 1 J„(s) m* + i Jkk(s) (4m%-3sf 


(67) 


and 


'q‘(«) =y{t-2mi){t- 2mi,) 4 (2 L) + LJ/S) - J 


3 1 


8 IGtt^ 


(1 + Lk) 


+ ml 


1 1 23 2 

32{-U, + LI--LI - Ll) + ^y- + 2 L^--L,) 


m„ 


16 LI + 32 


^9 9^"" 6 ^^ 6 ^" 2 ^" 


+ 1 S 


rjTT 


8 L^- -- 


1 1 


(1 + Lk) 


+ tm^ 


1 1 


m 


log(^) 


L3 167r2 
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B.2 The rjTT'^ K^K'^ amplitude 

The three functions involved in the NLO contributions to the amplitude were denoted 
as Uy, Wy and hhi. They can be expressed as 
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and 
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and, finally, 
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B.3 The amplitude K^K+ K^K+ 

The 0{p'^) contribntions to this amplitnde involve five fnnctions: 17i, Vq, Vi and 

Wq"^ (see eq. ( 11 )). Uq'^ can be expressed as 
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The function Ui{s) reads 
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Next, the functions Vo(t), Vi{t) read, 
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Finally, the function Wq^{u) reads, 
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C Asymptotic interpolation 

We describe here the simple ansatz which we used for interpolating the S'-matrix phases 
^ 11 , 622 and the inelasticity rj in the asymptotic region si < s < cx). Let F{s) be one 
of the functions (fii(s), ^ 22 ( 5 ) or arccos(? 7 (s)). We introduce a point S 2 = Si + e close 
to Si and we assume that F{si), F{s 2 ) are given. We denote S 3 = cx) and we use the 
asymptotic conditions (see (32)) 

(5ii(s 3) = 27r, 522(^3) = 0, r ]{ s 3 ) = 1 . (77) 
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Thus -^(ss) is also known. We introduce a function u{s) 


u{s) 


1 

1 + log — 

Si 


(78) 


which maps the range [si,cx3) onto the hnite range (0,1] and dehne F{s) through a 
simple Lagrange polynomial interpolation i.e. 




(u(s) 
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Uj){u{s) - Uk) 


(79) 


with Ui = u{si) and i, j, A; is a cyclic permutation of 1, 2, 3. 
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